Quantum critical behavior near a density-wave instability in an isotropic Fermi liquid 
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We study the quantum critical behavior in an isotropic Fermi liquid in the vicinity of a zero- 
temperature density-wave transition at a finite wave vector q c . We show that, near the transition, 
the Landau damping of the soft bosonic mode yields a crossover in the fermionic self-energy from 
E(fc,w) « E(fc) to E(fc,u>) « E(ti;), where k and ui are momentum and frequency. Because of this 
self-generated locality, the fermionic effective mass diverges right at the quantum critical point, not 
before, i.e., the Fermi liquid survives up to the critical point. 

PACS numbers: 71.10.Ay 71.10.Ca 71.10.Hf 71.18. +y 
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Introduction. An isotropic Fermi liquid may experi- 
ence various intrinsic quantum instabilities. They are 
characterized by divergence of the corresponding static 
susceptibility and emergence of a gapless bosonic mode 
in the collective, two-particle excitation spectrum. The 
instability point at zero temperature, occurring at a par- 
ticular value of a control parameter such as electron con- 
centration, is called the quantum critical point (QCP). 
Near QCP, the interaction between the soft bosonic mode 
and low-energy fermions often leads to singular behavior 
of the fermionic self-energy H(k,u>) and divergence of the 
fermionic effective mass m*. A well-known example is 
the divergence of to* near ferromagnetic instability 0, 
which occurs at the wave vector q — 0. 

In this paper, we study the divergence of to* in an 
isotropic Fermi liquid near a zero-temperature charge- or 
spin-density- wave instability occurring at a nonzero wave 
vector q c < 2&f, where k F is the Fermi momentum. We 
argue that, for q c ^ 0, the behavior of to* near QCP is 
rather tricky, and the analysis requires extra care. 

We consider a model in which fermions ipk interact by 
exchanging a soft bosonic mode V(q) 



Hint = - y 



q,k,k' 



(1) 



The soft bosonic mode is peaked at the wave vector q c 
and can be in either spin or charge channel 



V(q) 



t- 2 + (\q\-qc) 2 



(2) 



Here g is an effective interaction constant, and £ is the 
correlation length, which diverges at QCP as a function 
of electron concentration or pressure. 

Interaction between soft bosonic modes © was stud- 
ied by Brazovskii [2| in the context of a crystallization 
transition in an isotropic liquid. Dyugaev Q applied 
the model JD and (0 to explain enhancement of the 
effective mass and specific heat in liquid 3 He, arguing 
that 3 He is close to a spin-density-wave transition. Ref. 



utilized the Brazovskii model to describe a magnetic 
transition in MnSi, where a finite q c is likely caused by 
the Dzyaloshinskii-Moriya interaction. The model Q 
and |PJ also applies to itinerant electrons in the vicin- 
ity of a ferromagnetic instability - a small but finite q c 
appears there as a result of an effective long-range inter- 
action between fermions due to the 2kp Kohn anomaly 
@. Refs. |S Si 

proposed that the model QJ and (J2J 
can explain the enhancement and possible divergence of 
the effective mass observed experimentally in the two- 
dimensional electron gas (2DEG) JSL as we ^ as m 3 Hc 
films I n the scenario of Refs. USBj the instability 
at q c develops as a precursor to the Wigner crystal in 
2DEG or to the crystallization transition in 3 He films. 

The effective mass m* is extracted from the fermionic 
self-energy S(fe, u>) defined by the Dyson equation 
G^ 1 (k,uj) — iuj — £fe — where G(k,oj) is the 

fermion Green's function, and Sk — vp(k — fc^) is the 
bare fermion dispersion counted from the chemical po- 
tential. The derivatives of Ti(k,co) determine the renor- 
malization factor Z~ x = 1 + id^T, and the effective mass 
m*/k F — l/v* F = Z~ 1 (vF+dkT,)~ 1 , where vf and v* F are 
the bare and renormalized Fermi velocities. Divergence 
of m* can be caused either by id^Y, — > oo (and, hence, 
Z — > 0), or by — > — v F . In the former case, m* would 
diverge at QCP, but not earlier, while in the latter case, 
the divergence of to* would generally occur at a finite 
distance from QCP. There exist other scenarios 01 f° r 
the divergence of to* , which do not invoke a density- wave 
transition, but we will not discuss them here. 

The interplay between 9 W E and 9fcS depends on 
whether S(fc, w) predominantly depends on momentum 
k or on frequency uj. The two alternative scenarios for 
the model of Eqs. (JIJ and |J2J where E depends only on 
k or only on u were advocated in Refs. |(| 0, E] and 
Refs. (HQ, correspondingly. In this paper, we show that 
the behavior of S in an isotropic Fermi liquid near a 
density-wave transition is actually rather involved. At 
some distance from QCP, £(fc,o>) w However, the 
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FIG. 1: The solid circular line represents the Fermi surface. 
The ring of radius q c and width £ represents the effective 
interaction (Sj via a soft bosonic mode. A fermion with the 
momentum k close to /cf strongly interacts with the two "hot 
spots" obtained by intersection of the Fermi circle and the 
interaction ring. The vector components q± and gy are per- 
pendicular and parallel to the Fermi surface at the hot spots. 



frequency dependence of V, generated by the Landau 
damping, makes E(/c,a->) predominantly w-dependent in 
the immediate vicinity of QCP. We show that, when the 
fcrmion-boson interaction g is smaller than the Fermi en- 
ergy Ep ~ vpkp, the crossover from E(fc,u;) ss E(fc) 
to £(&, ui) ~ E(w) is separated from the weak-to-strong 
coupling crossover near QCP. The latter occurs when the 
dimensionless coupling A ~ (g / Ep)(^kp) oc £ becomes of 
the order of 1. On the other hand, the crossover from 
E(fc) to E(w) occurs at a small A — (g/Ep) 1 / 2 < 1, 
where E is still small, and <9fcE does not reach —Vp. 
Once E(fc, u>) becomes E(w), only <9 W E matters, i.e., 
m* = (Zvp)^ 1 diverges with Z _1 at QCP, but not ear- 
lier. We present calculations in the 2D case, but the 
results are qualitatively valid also in the 3D case. 

Momentum- dependent self-energy E(fc) away from 
QCP. In the Hartree-Fock approximation, the exchange 
diagram with the effective interaction J2J) gives 



E(fc,w) 



dn d 2 q 



(2tt) 3 
d 2 q 

W) 



G(k + qM + Lu)V(q) 
2 n F (e k+q )V{q), (3) 



where np(e) is the Fermi distribution function. The in- 
tegration over q in Eq. © is restricted by the conditions 
that the vector k + q lies inside the Fermi circle and the 
vector q belongs to the ring of radius q c and width £ _1 
centered at the vector k, as shown in Fig. ^ Clearly, 
E in Eq. J2J does not depend on u>, but it does depend 
of fc, because the area of the ring inside the Fermi circle 
changes with k. 

The derivative of Eq. © with respect to k, taken at 



k = kp, is given by the integral along the Fermi circle 
<9E f d 2 q 



dk 



■S(e k+q )v k+q V(q). 



(4) 



where 5(e) is the Dirac delta-function, and v k + q is the 
Fermi velocity at k+q. For large £, the integral l@J comes 
from the vicinity of the two "hot spots" q c obtained by 
intersection of the Fermi circle and the circle of radius q c 
centered at the point kp on the Fermi circle (see Fig.^l. 
Decomposing the deviation from the hot spot q = q q c 
into (q±,q\\), as shown in Fig.^ and integrating over 
first, we obtain <9fcE = fe<9fcE, where 

<9E m* v F 1 

— - = -A Vp cos 6>o, — = — = - r t-. (5) 

ok m vp 1 — A cos 0$ 

Here 9o is the angle between kp and kp + q c in Fig. ^ 
such that sin(6* /2) = q c /(2kp), and 



A 
2 



1 

Vp 



(2tt) s 



V(q) = 



9t 



Attvf cos(6» /2)' 



(6) 



We see that, if Acos#o > (which implies q c < ^/2pF for 
A > 0), the effective mass increases with A and nominally 
diverges at Acos^o = 1, while f is still finite, and QCP 
is not reached yet. 

Crossover to the frequency- dependent self-energy E(w). 
To verify whether Eq. © holds up to A ~ 1, where m* 
diverges, we need to go beyond the Hartree-Fock approx- 
imation and include the full fermionic and bosonic prop- 
agators and vertex corrections into the self-energy dia- 
gram. We assume and then verify that, for g/ Ep <C 1, 
the higher-order corrections predominantly renormalize 
V in Eq. while vertex corrections and the renor- 
malization of G can be neglected for arbitrary A. The 
renormalization of V(q) originates from the electron po- 
larizability 



H(q,0)= J 



d 2 kdn i i 

(27r) 3 i(uo + fi) - e k+q iuj - e k 



(7) 



via the relation V^ 1 (q, 51) = V^ 1 (q)+H(q, ft). For q near 
q c , the static part of n(g, f2) comes from the fermions 
with high energies and we assume that it is already in- 
cluded into Eq. |(2J), which implies that £ is the exact 
(renormalized) correlation length. The dynamical part of 
H(q c , ^) comes from low energies and describes the Lan- 
dau damping of the bosonic mode due to its decay into 
particle-hole pairs. For f2 <§C vpq c , lmll(q c ,il) oc 
Inserting the Landau damping term into Eq. @ , we find 



V(q,il) 



S- 2 + (q-q c ) 2 +j\ny 



7 



9_ 



(8) 



Re-evaluating E(fc,w) in Eq. J2J for the full V(q, il), we 
find that it now depends on both k and u>. Notice that 
causal analytical properties require that the interaction 
constant g in Eq. JSJ must be positive: g > |l3j . 
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We present the results for H(k,u>) first and discuss the 
details of calculations later. For small €k = vp(k — kp) 
and ui, we then obtain 

E(fc,0) = -A cos e e k h k (v), (9) 
T,(kp,uj) = -\iu)huj(rj), (10) 
V = jE F e~ \ 2 (E F /g). (11) 

Here and below we subtracted the renormalization of the 
chemical potential from E(fc, cj), i.e., redefined E(fc,w) = 
E(fc,w) — E(fcp,0). The functions hk(r)) and h u {r]) have 
the following asymptotic behavior. For 77 <C 1, frfe(?7) = 
1 + 0(rj) and hu,{rj) ~ 77 ln(l / 77) , i.e., the momentum- 
dependent piece in E almost coincides with Eq. (JjJJ, while 
the frequency dependence of E is weak. This is natural, 
because small 77 corresponds to small bosonic damping 
7. However, for r\ 3> 1, we find the opposite behavior: 
h k (rj) oc 77- 1 / 2 < 1 and h u {q) = 1 + O^ 1 / 2 ). In this 
case, the momentum dependence of E is weak compared 
with the Hartree-Fock approximation, while its frequency 
dependence is strong. Moreover, E does not depend on 
7 explicitly. Thus, the limiting forms of E(fc,w) are 



E(fc,w) 



-Acos6» £fc, 77 < 1, 
— Xiuj, 77 » 1. 



(12) 



When the system approaches QCP, and £ increases, 
the parameter 7/ cx £ 2 changes from 77 1 to 7/ ^> 1. The 
crossover between the two asymptotic limits in Eq. 112(1 
takes place at 77 ~ 1, which corresponds to 



A ~ A cr = Vs/^f < 1. 



(13) 



Thus, the upper line in Eq. 1(12(1 stops being applicable 
already at A < 1, before A can generate a divergence in 
Eq. ©. In the vicinity of QCP, the lower line in Eq. ltT2l 
applies, and, instead of Eq. Ip|l. we find 



m 
m 



1 

Z 



1 + A. 



(14) 



We see therefore that the effective mass in Eq. ((1411 di- 
verges only at QCP, where £ — * 00, but not before, con- 
trary to the conclusion one could draw from the Hartree- 
Fock approximation. This is the central result of the pa- 
per. Notice that the requirement g > 0, mentioned after 
Eq. ©, guarantees that Z < 1, because sgn(A) = sgn(<7). 

Further, since A cr <C 1, vertex corrections and renor- 
malization of the fermionic G in Eq. are small at 
A ~ A cr and can be safely neglected. This justifies our 
approximation of including only the renormalization of 
the bosonic propagator. Moreover this approximation 
actually remains valid even at larger A > 1. Indeed, the 
modifications to Eq. due to vertex corrections and 
residual momentum dependence of E are small in the pa- 
rameter y^g/Ep <c 1 and can be safely neglected even 
when A = 0(1). Although the fermionic E(w) is not small 
at A = O(l), using the renormalized Green's function 



G^ 1 (k,u>) = iu>(l + A) — VF(k — kp) = Z[iuj — v F (k — kp)] 
in Eq. (0 does not modify the lower line in Eq. i|12[l . 
because the extra factor Z and the renormalization of 
v* F = kp/m* compensate each other Similarly, the 
coefficient 7 in Eq. © does not change, because the fac- 
tor Z 2 coming from the two Green's function in the po- 
larization bubble Q compensates the renormalization of 
the factor l/vp — ► l/(v F ) 2 in the expression JHJ for 7. 
This behavior is typical for the Migdal-Eliashberg-type 
theories [TEEl- 

Anomaly in the calculation of self-energy. Now we 
present details of the self-energy calculation and also ex- 
plain why A cr vanishes if the fermionic bandwidth (~ Ep) 
is set to infinity. Linearizing the fermionic dispersion 
near the two hot spots, we introduce Q — vpq± and 
Ik = (fe — kp) ■ Vk F + Qa = efccos6*o, where the vec- 
tor kp is selected parallel to fc (see Fig. [TJ. Then 
E(fc, lu) = E(fe, lu) - E(fc F , 0) is 



J(fc,w) 



(iuj - e k ) I(k,u>), 

dnd( v(n,Q 
2tt [i(oj + n) - h - c] (*n-C) 



(15) 
, (16) 



where we introduced V(£l, £) similarly to Eq. (JSJ) 



Vp 



dq\\ 

W) 



X 



(17) 

Notice that, to this accuracy, V(Cl, C) does not depend 
on C, i.e., V(n,C) = V(Q). 

The evaluation of I(k,u>) in the limit k — > kp and 
lu —> requires care, because the integrand in Eq. 1(16(1 
contains two closely located poles separated by u> and Ik- 
If we approximate V(fl) by a constant V"(0) = A, then, 
nominally, the integral i(16|) is ultraviolet-divergent and 
depends on the order of integration over 11 and £. 

To evaluate the integral correctly, one must keep in 
mind that Eq. ((16(1 is approximate, and higher-order 
terms in {q—q c ) in G and V always make the integral over 
q convergent at q— q c ~ kp. If 7 = in Eq. ((17(1 . then the 
integral over O must be taken first, because its conver- 
gence is provided only by the fermion Green's functions in 
Eq. (fl?)|) . In this case, we obtain I(k,u) = \ek/(iuJ — e&) 
and E(fe, u>) — E^ = — Aefc, reproducing the top line of 
Eq. C2J. 

On other hand, if 7 is large in Eq. ((17(1 . then V(SY) 
strongly depends on fl and provides convergence of the 
integral over f2. In this case, it is appropriate to inte- 
grate over £ first, over the region where the linearized 
expression l(16|) is valid. Taking the integral over £ first, 
we obtain I(k, iuj) = Xiio/(iuj — efc) and E(fc, to) — —Xiu, 
reproducing the bottom line of Eq. ((12(1 . Notice that, al- 
though the frequency dependence of V(fi) is essential to 
determine the correct order of integrations, the strength 
7 of this dependence drops out from the final answer. 
This situation bears mathematical similarity to the chi- 
ral anomaly in quantum field theory |l7l [l8| ■ 
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The crossover between these two cases takes place 
when the characteristic ft in Eq. I|17|) becomes of the 
order of £ ~ Ep. Using the definition 1)110. we find that 
the cases of weak and strong frequency dependence cor- 
respond to 77 < 1 and 77 > 1, as in Eq. (|T2|l . 

The fact that the crossover occurs at n ~ 1, i.e., at 
A cr -C 1, is a consequence of V(q, Q) being peaked on 
a circle \q\ — q c . If V(q,SY) were peaked at a given 
vector q c in a crystal, then V would have the conven- 
tional, Ornstein-Zernike form V(Q, q~±) = f dq~\\/(£,~ 2 + 

3± + 9f + 7l^l) « (£~ 2 + <?! + 7l^l)~ 1/2 - In this case, 
the crossover takes place when all three terms become 
comparable: ~ q\ ~ 7|fi|. Since typical Q ~ vpq~±, 
the crossover in the vector case occurs at 7«f£ ~ 1, i.e., 
A C r ~ 1 C3- However, m* does not diverge even when 
E remains E(&) up to A ~ 1, because the correction to 
velocity <9fcE = — Xvk+q c is not antiparallel to vy, in the 
absence of nesting, so the magnitude of the Fermi veloc- 
ity does not vanish for any finite A [lij |. 

For completeness, it is instructive to see how the 
crossover from the top to the bottom line in Eq. I|12|) hap- 
pens if we always integrate over first in Eq. (|16|l . Let 
us deform the contour of integration over to either up- 
per or lower complex half-plane. For — ey. < ( < 0, each 
half-plane contains just one pole. Wrapping the contour 
around the pole and integrating over £ within the speci- 
fied limits, we obtain the E^ contribution to £. If V does 
not depend on Q, the calculation stops here. However, 
when V(tl) depends on Q and is given by Eq. fT7|) . we 
also need to consider a contribution from the branch cut 
in V(Q) along the imaginary axis of f2 where Q = iv + 5 
and |fi| = ivsgnS. Evaluating the contribution from the 
branch cut and combining it with the contribution from 
the pole, we find £ = — Ae^ + (iuj — ej)4 ci where 

/ bc - (2/tt) / d( / d V lmV{iv)/{v + Q 2 

Jo Jo 

2 P , f°° dvo ( A \ , x 

= -/ dz - Im -== . (18) 

K Jo Jo {Vj + Zy \y/l-ivoJ 

Here we introduced dimensionless variables z — C7^ 2 
and w = v^i^ 2 . For r\ <C 1, Eq. (|18|l gives a small 
/be ~ \r\\a.{\jr\). In the opposite limit 77 3> 1, the in- 
tegral over z can be extended to infinity, and the integral 
over xd yields /b c = V(Q) = A via the Kramers-Kronig 
relation. Notice that this result does not depend explic- 
itly on the detailed form of the frequency dependence in 
V(Cl) as long as the integral (|18(l quickly converges and 
can be extended to infinity. Substituting /b c into E, we 
reproduce both lines in Eq. |Q for 77 -C 1 and 7? 3> 1. 

Conclusions. In this paper, we studied the quantum 
critical behavior of an isotropic system of fermions near 
a T = transition into a density-wave state with a fi- 
nite momentum q c . We demonstrated that, upon ap- 
proaching QCP, the fermionic self-energy crosses over 
from S(fc, lu) w E(fc) to S(fc, w) w E(w). We showed that 



the crossover occurs while the dimensionless coupling A 
(which diverges at QCP) is still small. We found that 
the effective mass remains finite and positive away from 
QCP, and diverges only at QCP as m* cx 1/Z cx 1 + A. 

Our results apply to both charge- and spin-density- 
wave instabilities. In the latter case, spin-orbital interac- 
tion generally induces anisotropy in the spin space, e.g., 
for easy axis, the interaction in Eq. is mediated by z 
component of spins. This only affects the numerical co- 
efficient in E, proportional to the number of fluctuating 
spin components, but does not change the conclusions. 

We thank D. Maslov for useful discussions. The work 
was supported by the NSF Grant DMR-0240238 (AVC); 
by US-ONR, LPS, and DARPA (VMG); and by the NSF 
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